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SUMMARY 


The  objective  of  this  program  is  to  conduct  theoretical  and  experi¬ 
mental  research  to  determine  the  electromagnetic  scattering  from  hete¬ 
rogeneous  dielectric  bodies  as  individual  bodies  and  as  a  cluster  of 
bodies.  The  discrepancies  in  the  literature  regarding  the  singularity 
of  the  electric  dyadic  Green* s  functions  were  resolved  and  a  unified 
and  consistent  view  is  presented.  Compact  range  scattering  measurements 
at  1  GHz  were  successfully  performed  to  obtain  measured  data  to  validate 
the  numerical  analyses.  Extensive  computations  were  made  for  a  variety 
of  dielectric  scatterers,  including  a  one-foot  bird  at  1  GHz.  The 
agreements  between  measurement  and  computation  were  good  except  for  the 
resonant  sphere,  for  which  the  calculated  resonant  frequencies  were 
shifted  by  about  20  percent.  Various  numerical  techniques  were  investigated 
successfully  for  implementation  in  the  volume  methods  to  treat  symmetrical 
scatterers  through  use  of  symmetric  matrices,  the  use  of  banded  matrices, 
and  virtual  memory.  A  new  local-file  manipulation  technique  for  handling 
large  matrices  has  been  explored  and  found  to  be  potentially  useful. 

There  has  been  very  little  research  into  the  problem  of  scattering 
by  dielectric  objects  of  complex  permittivity.  Future  research  in  this 
area  should  include  the  investigation  of  the  surface  integral  equation 
technique  and  the  exact  solution  for  the  dielectric  prolate  spheroid. 
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INTRODUCTION 


During  the  past  decade,  there  has  been  a  rapidly  rising  interest  in 
electromagnetic  problems  involving  dielectric  objects.  The  interest  in 
this  area  arises  from  a  multitude  of  military  and  civilian  needs.  Human 
beings  are  increasingly  more  exposed  to  microwave  radiation  hazards  on 
aircraft,  ships,  military  installations  and  even  in  their  homes.  The 
knowledge  of  the  electromagnetic  scattering  from  birds,  animals  and 
humans  is  essential  in  the  analysis  of  detection,  identification  and 
interference  problems  in  radar  systems.  Biological  and  medical 
appl ica t ions ,  such  as  blood  thawing,  enzyme  inactivation  and 
hyperthermia  treatment  of  cancerous  tissues,  also  demand  accurate 
knowledge  of  the  electromagnetic  fields  in  dielectric  bodies.  Thus, 
e 1 ec tromagnet ic  scattering  by  dielectric  bodies  is  a  fundamental  and 
important  problem,  and  solutions  for  these  scatterers  have  wide 
appl i cat  ion . 

Since  September  1978,  Georgia  Tech  has  been  supported  by  the  Deputy 
for  Electronic  Technology  (RADC/EEC),  Air  Force  Systems  Command,  under 
contract  FI 9628-78-0-0223  to  conduct  a  two-year  research  program  in  this 
area.  The  emphasis  is  focused  on  analyses  and  measurements  of  the 
scattering  characteristics  of  heterogeneous  dielectric  objects  as 
individual  bodies  and  as  a  cluster  of  bodies.  Progress  to  date  has 
included  the  clarification  of  discrepancies  in  the  literature  on  the 
singularity  of  the  Green's  function  in  the  source  region,  and 
measurements  and  calculations  of  the  scattering  cross-section  of 
dielectric  bodies  of  arbitrary  shapes  and  complex  dielectric  constants. 
In  addition,  various  numerical  techniques  have  been  investigated  and  a 
new  local-file  manipulation  technique  was  explored  and  found  to  be 
po  t  e. r.  t  i  a  1  1  v  useful  . 

Numerical  techniques  developed  in  the  current  research  program  are 
capable  of  producing  fairly  accurate  data  for  objects  less  than  one 
free-space  wavelength  long.  There  are  models  for  which  the  present 
technique  is  highlv  accurate,  and  there  are  geometries,  such  as  the 
sphere,  tor  which  the  present  technique  is  not  quite  satisfactory 
probably  due  to  the  inherent  deficiencies  of  the  volume  integral 


The  success  of  compact-range  scattering  measurements  at  a  frequency 
of  1  GHz  represents  an  advance  in  the  state-of-the-art  of  scattering 
measurements  using  the  compact  range  technique.  This  extension  of  the 
compact  range  scattering  measurements  to  lower  frequencies  demonstrates 
the  versatility  and  usefulness  of  the  compact  range  as  a  general  purpose 
EM  measurements  facility. 
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I.n  the  source  region  where  r  -  r  ’  ,  questions  arise  as  to  whether  1-4)  is 
stiLl  valid  and,  if  valid,  how  to  compute  it  sinev  tiie  i<  ’  term  in  (a  ! 
is  in  general  not  integrable  jjj.  For  the  static  case  in  which  k  -  u, 
the  rigorous  treatment  by  Kellogg  [3]  is  applicable.  Using  the  method 
of  retarded  potential.  Van  Biadel  (4 \  showed  that 

f  d  ( r ) 

K(r)  PV  I  J(rr)  •  i.;1'  (r,rf)  dV' - .  (h) 

J  V  "  ”  .  3j  • 


”PV"in  Lquation  (b)  denotes  ’’Principal  Value”  integration,  which  specifies 
an  integration  carried  out  over  the  volume  V-V  ,  where  V  is  an  infinitesimal 
spherical  volume  centered  at  r  . 

With  a  rigorous  classical  procedure  similar  to  Kellogg’s  treatment 
of  the  static  case,  Fikioris  [5]  derived  the  following  expression 
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where  a  is  the  radius  of  a  finite  spherical  volume  centered  at  r. 

He  also  indicated  that  does  not  have  to  be  spherical  but  is  not 
permitted,  in  general,  to  be  infinitesimal.  If  assumes  other 
geometries,  the  third  term  on  the  right  side  of  Equation  (7)  will 
have  to  be  modified  accordingly.  Equation  (6)  also  answered  implicitly 
the  question  as  to  how  to  evaluate  the  principal-value  integration  in 
Equation  (6).  Thus,  for  a  spherical  principal  volume  of  radius  a  and 
a  constant  current  J, 


J (r)  •  G  (r , r  ' )  dVf 


e  ^a(J  +  jka).  (8) 


As  will  be  discussed  later.  Equations  (6)  and  (8)  form  the  basis  on  which 
numerical  analyses  using  free-space  electric  dyadic  Green’s  function 
were  carried  out  [ 6—8 j . 

Chen  [9]  showed  that  the  infinitesimal  principal  volume  V  can  assume 

other  geometries  as  long  as  the  integrals  over  V  and  V-V  are  properly 

c  c 

handled.  He  also  derived,  as  examples,  expressions  similar  to  Equation  (6) 
for  principal  volumes  in  the  shapes  of  the  cube  and  the  circular  cylinder. 

For  fields  in  the  bounded  regions,  the  electric  dyadic  Green’s  functions 
for  rectangular  cavities  and  waveguides  were  discussed  by  various  authors 
1 10-16 j.  Their  expressions  are  in  the  form  of  eigenf unction  expansions 
plus  a  term  containing  the  Dirac-delta  function.  Perhaps  considering  the 
deita  function  to  be  the  only  singular  term,  Yaghjian  [ 1 7  J  explained  the 
difference  in  the  delta-function  terms  between  Tai,  et  al .  [11 j  and 
Rahmat-Samii  [12]  as  being  due  to  their  different  choices  of  the  principal 
volume  V  and  emphasized  the  need  to  include  in  G  the  shape  of  the  principal 
volume  involved.  Johnson,  et  al.,  [16]  also  pointed  out  the  ambiguities  in 
the  literature.  Recent  l  v,  Lee,  et:  al.  1 18 1  examined  this  problem 
1 rom  the  viewpoint  of  the  potential  theory  and  obtained  expressions  for 
the  genera i  problem.  Most  of  their  conclusions  support  existing  results. 


buL  questions  concerning  the  usefulness  of  an  ini initesimal  principal 
volume  V  in  numerical  analyses  were  raised. 

Proper  handling  of  the  electric  field  integral  equation  in  the 
source  region  is  essential  when  using  it  in  numerical  analyses  involving 
dielectric  scatterers.  The  difficulty  arises  in  the  calculation  of  the 
"self-cell*1  or  self-coupling  matrix  elements  that  must  be  generated  when 
using  the  method  of  moments.  For  the  free  space  [ 6, 8]  and  half -space 
[7],  Equations  (6)  and  (8),  which  are  consistent  with  Equation  (7),  were 
employed  in  the  scattering  analysis  of  arbitrarily-shaped  dielectric 
bodies.  For  the  rectangular  waveguide.  [19],  the  expressions  of  Tai  [10] 
and  Rahma  t-Sami  i  [12]  were  used.  To  the  practicing  engineer,  it  is 
desirable  to  remove  all  discrepancies  and  subtleties.  in  the  following, 
it  is  shown  that  uniformity  and  consistency  in  this  subject  can  be  reached 
*n  the  basis  of  mataemat leal  analyses  and  numerical  experiments. 

ij_.  A  *.  a  !_f  j  .si  and  >_aa  ;  stun t View 

Discussions  of  me  singularity  associated  with  the  electric  dyadic 
iti'i-on's  :  u.u;  t  i  .*n  in  List  free  spruce  are  comprehensive  in  the  literature 
mu  :.ons  are  not  a. rious.lv  divided.  However,  t  .?j  th>j  bounced  regions, 
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rt.tan^uiar  v.ivt?  t'u  1  du  a  no  c«i  v  i  l  v  ,  and  the  hair  splice  ,  t  t  ,  nas  i'lcii 
we  1  *  documented  in  the  literature.  Since  the  ini  »*.•:  in.;  :  fmiat  ■;  mi  <  -c  i 
is  iess  singular  than  K  ,  in  ract  it  is  K  in  tree-space,  tne  .inttmrai 

in  Kquation  (9  b)  converges  13,  p.Mfcj.  bv  m::  a  1 1 1  ut  .  :m  count  i  i  *  n  •  co  .• 
into  he  ua  t  ion  l9a)  and  as  sum  me  that  .J  satisfies  the  so-called  hoider 
Condition",  it  was  proved  bv  Kellogg  (.3}  and  Kikioris  15)  for  k  =  0  and 
k  *  0  respectively  that  VV  •  A  exists  for  the  tree-snace  case.  Tins 
conclusion  was  later  expanded  to  the  bounded  regions  of  waveguides  and 
cavities  by  Yagjhian  1 1 7  J  and  Lee,  et  al.  [18].  It  is  interesting  to  nut 
that  by  assuming  the  first-order  differentiability  of  J,  Van  Bladel  [  a  j 
was  able  to  establish  similar  results  bv  usine  the  relations  between  suri 
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and  volume  integrals  such  as  the  divergence  (Gauss)  theorem.  The  Holder 
condition  is  stronger  than  continuity,  but  may  be  either  weaker  or  stronger 
than  differentiability. 

The  potential  approach  is  based  on  familiar  and  well-established 
classical  analytical  techniques  and  there  are  no  identifiable  disputes 
regarding  it.  If  we  choose  the  delta-function  approach  based  on  Equation 
(3),  we  should  treat  it  as  a  distribution.  gP  is  no  longer  a  classical 
function,  but  is  defined  only  when  it  is  used  in  an  integral  [2,21].  Now,  if 
Gf  is  expressed  in  terms  of  functions  that  appear  to  be  conventional,  one 
may  be  misled  to  regard  ge  as  being  also  conventional.  For  example,  [22,  p.  20CT 


2 


(jVd-  }  +  k“ 

ti  ?T  K 


the  left  side  of  the  equation  appears  to  be  a  classical  function  but 
can  only  be  defined  as  a  distribution  at  R  -  0.  Note  that  the  left  side 
of  Equation  (11)  contains  a  R  singular  terra  and  is  closely  related  to  the 
singular  behavior  of  the  free-space  electric  dyadic  Greenfs  function. 

his tor icail} ,  the  tree-space  G  is  closely  associated  with  the 
method  of  potential  and  G  for  the  bounded  region  is  more  closely 
associated  with  the  method  of  distribution.  They  are  discussed 
sCf>~*raLei'  as  iol.1  o  w  s  . 
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i*  tae  current  distribution  satisfies  the  so-called  Holder  Condition, 
the  eiectiiv,  in  the  tree-space  can  be  expressed  in  terms  of  the 

eieeir i ivauic  green’s  function  in  the  following  general  form 
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on  trie  Hie t hoc  oi  potential,  using  relations  between  surbice  and  volume 

integrals,  suen  as  the  divergence  theorem  [2oi.  'inis  mternutivi  me  mod  is 
similar  to  that  of  Van  bladei  [4)  and  the  classical  portion  oi  Let-,  et  ax. 
f lb j  and  must  rely  on  the  existence  oi  the  iirst  derivative  oi  the  current 
J  at  the  singular  point  r  =  r’.  In  practice,  this  slight  ditiereneo  m 
requirements  tor  J  is  not  considered  to  be  significant. 

If  we  let  the  maximum  cord  of  V  Vanish,  the  following  expression 
by  Yaghjian  [ 1 7  j  emerges  from  Equation  (12) 


* 


as  well.  It  must  also  be  pointed  out  that  the  evaluation  of  the  integrals 
over  V-V^  and  V£  for  a  finite  is  by  no  means  simple  and  is  likely  to 

depend  on  the  particular  problem  under  consideration.  Thus,  the  claim 
that  an  infinitesimal  principal  volume  is  not  suitable  for  numerical 
computation  appears  at  best  premature. 

e 

b .  Q  in  Regions  Bounded  by  Conducting  surfaces 

0 

A  very  important  feature  of  G  for  the  bounded  region  is  that  ex¬ 
pressions  in  the  literature  have  been  derived  from  the  Ohm- Rayleigh  method 
which  expands  the  Dirac^delta  function  in  terms  of  the  eigenfunctions 
of  a  homogeneous  Helmholtz  equation  [2].  Consequently,  GC  derived  by 
this  method  is  based  on  the  distribution  theory  and  must  be  understood  as 
a  distribution,  or  generalized  function.  Furthermore,  being  a  distribution, 
Gt  does  not  necessarily  need  the  special  handling  for  its  singularity 
through  the  prineipai  volume  approach  since  the  purpose  of  the  distribution 
expressions  are  tor  the  handling  or  singular  integrands.  Both  the 
distribution  approach  and  the  principal  volume  approacn  are  discussed  in 
Luc  ioliowing. 

iagnjian  il/j  anti  Lee,  et  al,[.lr  (  eutlinec  precis  n«  snow  that  the 
singularity  associated  with  trie  electric  d.*io:c  wivtc’s  lunctien  in  .» 
bounced  resnn  is  exact  lv  the  same  as  that  for  the  tree-snare. 

Snecificni.lv,  i quat  i  uic  til*  and  ■!  >'  ••  }  :  :  : b  am  led  re.:  an¬ 
as  well  if  (•*"  S']  t  i  5  1  i  ex  the  ruiiarion  on::i  i  r  i  on  n  >:  .  w  :  ; 

is  either  ir.rir.ite  or  per:  ect  tv  concur'  :  n  r. ,  :  n  .  *uw  c  .  i •  ■-  an. 

trie  c  o :  i  d  * t  i  v  1 1  v  the.  b-aiuu  .rv  :  ^  ;  an?  :  •  t  ,  :  .  .-*k  sti.l  •:  an  1  <.  :  ■ 

ust'  the  equivalent  volume  current  t  necr  t  r:  t recast  tat  p  :*  ■  n  ,  ■  r.  :  i  a  :  :u 
case  involvin'*.:  either  in: ini  tv  renaa.a  ar  :n-r :  o  r  :  i  v  rvma  Lai. 


1  U  e  a 


.  }vl  *  . r. in:;  of  Equu.  '•  :  (12  :  •  or  the.  ui:2  ■  ..  ’.’ti-  in  .  ^rr  : ».  a 
o:.iL  nv  metnoa  of  r  i .. !  ..me  the  v..-  .  .  ■  pro  [  ■'  ■* 

; '  ;>je'  •  ‘  •  ■'  '  ‘  •■'•:•  ■■■'.. 

b  e  differentia  b  1 »;  . 

Since  G~  lot*  trie  bounded  regions  is  often  derived  from  the  0 m — Raleigh 
method,  and  hence  is  a  distribution,  it  can  he  used  either  in  conjunction 
witn  Equation  (12)  or  directly  as  a  distribution.  By  either  method,  one 
snouid  obtain  trie  same  results  for  trie  electric  field. 

Numerical  computations  involving  G  in  a  bounded  region  were  ca' inert J 
::■>  Wang  LiVj  in  determining  the  field  intensity  in  an  arbitrarily-shaped 
dielectric  or  Dioiogical  body  inside  a  rectangular  wa/eguide.  Since  G 
i derived  by  the  .ihni-Ka  Le  i  ch  method ,  trie  question  as  to  what  G  is  at 
r  -  r  :  need  not  he  a  mat  tel  .  :  ncerr; .  Me  a  awn  lie ,  trie  principal  tw. '  unc 

;  at  ••  n  ;  x  1 1  .ni  l  n.nn  ;ue  is  r.n  .  i  n .  *  e  a:.;  use:  u  L  c  ;  trier  as  a  principal  or 

,  i  ary  too  i  out  is  no  .oncer’  as  1 1  m  i  or  the  3  ret  space  regie 

me  most  severe  a: i ticu.lv  encountered  in  the  anal vs  is  ot  the  waveguid 
pro h gen.  in  reieience  i^  is  the  nanJiiiig  oi  trie  singular  integrand,  w’hich 
was  only  hi  iei  iv  uescriheu  in  tne  paper.  ine  method  was  supported  bv 
o>  tn  L  he  uist  r  .  i)uL  Lea  meorr  aau  the  principal  volume  approach.  In  the 
i  liter  meinod,  a  ‘piii  pox'  principal  was  cnosen  but  other  geometries 
couiu  oe  used  as  long  as  trie  integrals  were  correctly  handled. 

me  electric  dyadic  Green's  function  pertaining  to  a  rectangular 
waveguide  was  first  derived  by  Tai  [10]  as  follows 

Gr,. 

where  u  is  a  double  infinite  trigonometric  series  resulting  from  an 
e i centime t ion  expansion  of  the  boundary-value  problem.  In  using  the 
principal  volume  metiioc.  w-  choose  a  pill  box  to  enclose  the  singular 


point  r_  as  shown  in  Figure  2.  Based  on  the  results  of  Yaghjian  [17],  we 
have 
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J  Z 
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(18) 


The  integral  over  Z1  in  Equation  (18)  can  be  integrated  and  is  described 
by  I  in  Equations  (19-20)  of  Reference  19.  The  integration  over  x  and  y 
leads  to  double  infinite  series  identical  to  those  in  Equations  (16-18) 
of  Reference  19.  The  series  did  not  appear  to  converge  when  it  was  first 
computed  numerically.  It  was  found  later  that  the  rapidly  oscillatory 
part  of  this  series  can  be  summed  up  in  a  closed  form  and  the  numerical 

value  of  the  resulting  expression  converges  satisfactorily  [19]. 

c 

We  can  also  view  G  as  a  distribution  since  it  was  derived  from  the 
delta-function  method  based  on  Equation  (3) .  Using  the  equality 


gn( x  +  t)  dt 


g  (x  +  t)  dt 
n 


(19) 


which  holds  for  every  convergent  series  of  distribution  [24,  p.  47], 
we  can  apply  it  successively  for  n  and  m  to  obtain,  for  a  constant  J, 

r  Z? 

H  H  £  f  U)  g  (y)  /  “  h(z,z’)  dz’ 

m  n  J  Zf 

J  Z 
z 

jeoc  (20) 

where  P  is  a  vector  determined  by  the  polarization  and  amplitude  of  J. 
Again,  the  integral  over  z!  and  f  and  can  be  summed  up  to  produce 
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lively  i:.q  uu  t.  re  ns  (  1 9-  ;  a!  u  in  !<*-•:  •_*  :viu*o  V.  Tliu.-; , 

(dm.  .  :  »/e*i  ..  r-  the  d:-  “Lijiir.ii.JU  *  ry  i  J  , 

v-  ■  :  cue  •.  \u,  r.>.h’!i  . 

cue  convergence  condition  required  .In  the  application  oj  .yum -on  (jwj 
According  to  Schwartz  [2T,26| 


V 


.  U  X. 


: ;  a;iu  cum 


~e  is  an  iu;-.-cer 


l  i  s  i  i  .*  o  cm  l  n  t  o  l  c  s  t  to  see  what  happens  i  £  we  choose  a  dir  i.  e  it  i;i. 

:t  i  :i..i  :m  i  ve  mum .  snort  as  a  i;ubo,  :  n  using  Kcjuatior.  i  1 7  )  .  in  t;iis  case. 

tne  see  one  tern;  on  tiie  right  oiue  ui  ric;  nation  t ib )  will  nave  to  be  changed 

to  -•  .J  *  ...  :  a*,  sordine  c  o  heuat i on  i.  L  3  )  .  .  t.  the  sum  time.  V  in 

r ioure  _  will  nave  to  be  reduced  and  the  integral  over  tne  added  volume 

in  v  must  oe  computed  and  added  to  tiie  integral  on  tin  loll  ot  1-quat.ion  (\ 8 )  . 

i 

wili  become  more  evident  in  the  next  section,  G*  is  singular  and  the 
integration  over  these  added  infinitesimal  volumes  does  not  vanish  and 
will,  in  combination  with  the  rest  of  the  terms,  yield  the  correct  vaLue 
for  the  electric  tieid. 

For  J  that  is  not  or  is  not  assumed  to  be  constant.,  liquation  (in)  can 
;>e  used  in  conjunction  with  liquation  (ib)  (with  G  replaced  dv  G  ,  oi  course) 
lor  the  principal,  volume  approach.  For  the  distribution  approach,  tin- 
technique  oi  equation  (ibj  can  also  be  used  in  Kquattoa  iT  >  .  i:n-  i  r:  >.*■■  era  l  j  on 
cl  Lin-  !  onutiUit  term  cor  respond  ing  to  i‘r>  has  noer.  discussed  arm  i.:u-  u  r;'i 
■ -i  resnonu  i t  t  ir.r’j  will  no  an  ordinary  initriai 
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Various  Forms  of  G  for  the  Rectangular  Cavities 

g 

As  has  been  pointed  out,  G  for  the  bounded  region  derived  by  the 
Ohm-Raleigh  method  is  a  distribution  and  must  be  so  treated.  Yaghjian  [17] 
and  to  some  extent  Johnson,  et  al.  [ 16  J  have  overemphasized  the  delta 
function  term.  In  this  section  it  will  be  shown  that  the  expressions 
of  Tai  and  Rozenfeld  [11]  and  Rahmat-Samii  [12]  are  mathematically 
identical  and  no  discrepancy  appears  to  exist.  This  finding  was 
previously  outlined  in  a  symposium  paper  [27]  and  the  full  derivation 
follows . 

The  electric  dyadic  Green’s  function  derived  by  Tai  and  Rozenfeld 
[11]  is  of  the  following  form: 
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The  expression  by  Rahmat-Samii  [12]  had  an  error  in  sign  resulting  from  a 
misprint.  After  correction,  his  expression  reads,  in  terms  of  the  e*^*’ 
convention  chosen, 
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Substitution  ot  Equations  (23-32)  into  Equation  (33)  yields 
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Substitution  of  Equations  (38-40)  into  Equation  (37)  yields  an 

0  0 

expression  of  £T  identical  to  Q.R  of  Equation  (33),  with  the  understanding 
that  ra  and  n  are  exchanged  and  that 
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if  at  least  two  of  the  three  integers,  m,  n  and  &,  are  nonzeros  (when  two 
or  more  of  them  are  zero,  the  series  term  is  zero).  We  can  also  prove 
directly  that  Equations  (22)  and  (33)  are  identical  either  by  substituting 
only  Equation  (40)  into  Equation  (37)  or  by  using  Equations  (38),  {39) 
and  the  following  relations 
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which  were  originally  employed  by  Tai  and  Rozenfeld  and  can  be  directly 
derived  from  well  known  identities  [20,  p.  581].  This  finding  that  the 
expressions  of  Rahmat-Samii  and  Tai  and  Rozenfeld  are  identical,  even 
though  they  appear  in  different  forms,  has  also  been  observed  by  Professor 
Tai  [21].  For  a  more  detailed  description  of  the  proof,  the  reader  is 
referred  to  an  earlier  interim  report  [28], 
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first  implemented  in  X  and  bands  [29]  »  and  was  re-eiu  lv  extended  to 
30  GHz  [30].  At  frequencies  below  2  GHz,  a  number  of  technical  dif¬ 
ficulties  ’nave  3 ong  been  general  1  y  recognized.  >L»  j,  r  y:  ■  La  .;es  in- 
• :  1 i 1 c  r.  •. : * ♦. .  e d  a d  i  :  f  r  a c  t  i  n n  o  f  •: h t*  r  e-  f  .1  c*  c  t  c;  r  ,  mu  1 1  i  pare  c  *  t  h e  i  j  ]  uni  n a  tine 
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Luned  source  with  an  output  pow’er  of  20  to  200  mw  was  built  for  the  measure¬ 
ment  system.  When  connected  with  a  lock-box,  a  frequency  stability  of  1 
part  in  10^  was  maintained.  The  dual-channei  phase  and  amplitude  receiving 
system  is  shown  in  Figure  >. 

The  compact  range  used  in  this  study  consists  of  a  12-foot  high  b\ 
lb- foot  wide  reflector  as  shown  in  Figure  h.  The  reflector  was  fed  bv  a 
24  in.  x  32  1/2  in.  rectangular  horn,  shown  in  Figure  /,  located  at  a  focal 
distance  or  about  12  feet.  Both  tne  reflector  and  horns  are  fabricated  by 
standard  methods  with  average  mechanical  tolerances, 
a  st  vrofoan  support  tor  t.h>-  target  under  test. 
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Figure  3.  Block  diagram  for  1  GHz  compact  range  RCS  measurement. 


Figure  4 .  The  arrange men 
scattering  men: 


At  1  GHz  and  lower,  the  multipath  propagation,  leakage  and  parasitic 
excitation  along  the  source,  cable  and  the  components  and  equipments  can 
cause  severe  difficulty.  In  fact,  leakage  from  the  source  was  initially 
found  to  be  only  30  dB  below  the  illuminating  field  in  the  quiet  zone. 

These  problems  were  overcome  through  careful  shielding  and  the  use  of 
microwave  absorbers. 

The  field  in  the  quiet  zone  was  probed  with  a  dipole,  and  amplitude 
variations  of  less  than  +1  dB  and  phase  variations  less  than  +5  degrees 
over  an  area  of  5  ft.  x  3  ft.  were  achieved.  Since  the  horn  used  is  an 
ordinary  rectangular  horn,  further  improvement  in  the  quiet  zone  illumi¬ 
nation  may  be  achievable  by  using  low-side-lobe  feeds  such  as  a  corrugated 
horn. 

The  ultimate  criterion  for  the  radar  cross-section  range  using  the 
cancellation  method  is  the  stability  and  depth  of  the  null  achievable. 

The  deeper  the  null,  the  smaller  echo  return  the  system  can  detect.  Also 
a  stable  null  insures  accurate  and  consistent  measurements.  To  achieve 
stability  and  depth  for  the  null,  sufficiently  high  power  and  frequency 
stability  of  the  source  are  essential.  The  sensitivity  of  the  receiver 
is  usually  sufficient  since  the  environmental  noise  in  the  range  is 
usually  quite  high.  In  the  1  GHz  measurements,  we  were  able  to  obtain 
a  null  depth  of  -50  dBSM  during  the  day  and  -60  dBSM  in  the  night,  which 
could  be  maintained  for  an  average  duration  of  1.5  to  2  minutes. 

The  sensitivity  of  the  compact  range  is  displayed  by  measurements 
shown  in  Table  I  on  small  conducting  spheres  shown  in  Figure  9  whose 
echo  areas  are  accurately  known.  The  close  agreement  shows  that  accurate 
measurements  can  be  made  for  small  scatterers  with  low  echo  return.  Figure 
10a  and  10b  show  measurements  on  a  circular  cylinder  2.76  wavelengths  long 
for  E  and  H  plane  aspect  angles.  Figure  11a  and  lib  show  the  measured 
data  for  flat  conducting  plates.  Figures  12a  and  12b  show  the  measured 
data  for  rectangular  conducting  boxes.  All  these  measurements  are  in 
good  agreement  with  data  in  the  literature,  as  can  be  seen  in  Figures  9 
through  12.  However,  there  is  some  confusion  in  the  literature  concerning  the 
polarization  of  the  data  which  remains  to  be  clarified. 
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'I  ABLE  I 


COMPARISON  BETWEEN  COMPACT  RANGE  MEASUREMENT 
AND  EXACT  CALCULATIONS  FOR  CONDUCTING  SPHERES 
CALIBRATED  WITH  THE  0.04^03*  SPHERE 


Radius  in 
Wave  1 eng  th  s 


RCS  in  dSSM 


Theore  t ical 


Measured 


04403 

-46.14 

-46.1 

05503 

<r 

c 

i 

-40.4 

06054 

-37.63 

-37.6 

06604 

-35.45 

-35.5 

07154 

-33.42 

-33.4 

07705 

-31.80 

-3:  .8 

104  56 

-25.44 

-25.4 

,11010 

-22.82 

-22.8 

Stainless-steel  spheres  used  In  the  calibration  of  the  scattering  measurement. 


Theoretical  (Oshiro) 


PLANE  ASPECT  ANGLE  a (DEGREES) 

-plane  pattern 

compact  range  measurements  and  other  ki 
Linder  (continued). 


Theoretical  (Oshiro) 


PLANE  ASPECT  ANGLE  a  (DEGREES) 


BACK  SCATTERING  CROSS  -  SECTION  (<JB) 


ELEVATION  ANGLE  (DEGREES) 

a)  Backscattering  pattern  as  a  function  of  elevation  angl 

Figure  11.  Comparisons  of  measured  and  calculated  backscattering 

cross-section  of  a  square  plate. 
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a/X 

b)  Backscattering  cross-section  as  a  f unction  of 
plate  size. 


Figure  1L. 


Comparisons  of  measured  and  calculated  backseat tering 
cross-section  of  a  square  plate. 


RADAR 


Figure  Id  .  Compact  range  scattering  measurement  oi  conducting  boxes 


3 i mu 1  a  t i o n  o_l  Biol i yg  i c a_^_  Hod i  <.-  * 

A  major  di:  l icuitv  in  measuremen mv*  ivinr  a  i. :  !  a:  -a]  ;s  t  no 

i  a ..  .  c onsistenc y  a n d  s ■  A  i  i  v  \  b i  r  d  or  an  i  na  :  unci r  a u e a 

priory  Logical  changes  a  i  -  tiu  -  ’.me.  Tny-*  may  .  •  a  less  today,  an  a  r  :.c  y 
may  move  during  r  he  tests.  inis  cii  I  :  •  c  ui  t  ■-  can  no  vt-rv  frustrating 
and  the  interested  reader  may  consult  reference  32  for  methods  considered 
by  Blacksmith  and  Mack.  [33]. 

Since  the  physical  configuration  and  the  dielectric  property  are 
the  only  parameters  involved  in  tne  scattering  measurements ,  it  is  feasible 
to  perform  the  test  on  a  simulated  model.  While  the  simulation  model 
removes  the  inconsistency  and  instability  inherent  in  real  biological 
bodies,  there  exists  some  degree  of  uncertainty  as  to  how  close  the  simu¬ 
lation  can  be  realistically  accomplished.  For  example,  it  is  difficult  to 
simulate  tne  feathers,  .skin,  bieoa  vessels,  etc.  However,  1  trainers  a  no  smn 
have  low  dielectric  constants  and  can  sometimes  be  conveniently  ignored, 
i’he  blood  vessels  have  htgn  water  content  but  are  usually  surrounded  musci 
t  is  i,  ue  s  wa  i  ch  a  iso  cental  r.  oat  o  r . 

Simulation  techniques  for  biological  bodies  were  developed  bv  uuy  |i^j, 
who  employed  various  chemicals  to  simulate  the  complex  dielectric  constant 
of  rat,  bone,  and  muscle  tissue.  His  method  was  used  in  the  fabrication  of 
bird  models  in  this  research  program. 

Before  making  the  model,  a  styrofoam  mold  is  constructed  to  hold  and 
supper*,  tne  moist  jellied  ,!Super-stu£  i  ”  plastic.  Trio  mold  is  derived  rrom 
a  biro  model  made  of  plaster  of  paris.  Figure  LJ  shows  the  side  and  i rent 
views  of  a  sitting  bird  model  made  of  piaster  of  paris.  Figure  I -  shows  try- 
side  and  front  views  of  a  flying  bird  model  made  of  plaster  of  paris.  Figure 
16  shows  the  side  view  ot  a  ”Super-stuf f ''  simulated  flying  Green  Wingtail,  SB 
in  a  st  v  ro : oam  b o 1 d e  r . 

The  ’'Super-stuff"  simulated  muscle  tissue  is  composed  of  the  following 
rna  L  e  r  i  a  1  s 

S.i  1  ine  solution  ( 1  2  g  sa  1 1  / 1  i  t  e  r )  76.3  ( b v  we  i  p  h  t  ) 

Powdered  p  o  1  y  e  t  h  y  j.  e  n  o  1^.2 " 

Sue f'”S  L  uf  f  3.4  7  . 


Fi^un 


(a)  Side  view 


(b)  Front  view 

L 3 .  Side  and  front  views  of  a  sitting  bird  made  of  plaster 
of  paris. 
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The  "Super-stuff”  is  a  jelling  agent  manufactured  by  Oil  Center  Research 
Corp.  in  Laffayette,  Louisiana.  The  mixing  process,  which  was  improved 
here  by  trial  and  error,  is  critical  to  the  homogeneity  of  the  simulated 
tissue.  Spectroscopic-grade  salt  is  first  added  to  deionized  water  in  a 
blending  mixer  in  an  oven.  After  reaching  200°F  the  solution  is  stirred 
for  about  2  minutes.  Fine  polyethylene  powder  is  then  slowly  poured  into 
the  solution  which  is  now  being  stirred  at  high  speed*  After  half  of  the 
polyethelene  powder  is  poured  in,  the  rest  of  the  powder  is  mixed  with  the 
super-stuff  and  poured  into  the  solution  being  stirred  at  high  speed.  The 
temperature  and  stirring  help  to  remove  bubbles  and  attain  homogeneity.  The 
temperature  is  then  raised  to  450°F  for  two  minutes  and  the  mixture  is  then 
allowed  to  cool. 

Although  it  is  usually  possible  to  make  a  simulation  model  to  meet  the 
required  dimensions,  it  is  not  easy  to  obtain  the  required  complex  dielectric 
constant  with  high  accuracy.  In  order  to  insure  reasonable  accuracy  in  the 
model,  the  in-vivo  probe  measurement  technique  [35]  was  used  to  determine  the 
complex  permittivity  of  the  model.  Ordinarilly  there  is  about  5  percent 
error  in  this  dielectric  constant  measurement.  In  ail  the  simulation  models 
measured,  the  disagreements  between  the  anticipated  and  measured  complex 
permittivity  were  mostly  within  5  percent  of  each  other. 


SECTION  [V 


An  exa:t  solution  for  turee-dimensionai  dielectric  scatterers  in  free 
space  exiscs  only  for  the  sphere.  For  scatterers  oi  arbitrary  siiapes, 
numerical  analyses  employing  the  volume  integral  equation  have  been  con¬ 
ducted  [6-7,36-38].  There  are  other  numerical  and  approximate  methods 
which  are  devoted  to  the  estimates  of  SAR  (Specific  Absorption  Rate)  [38], 
which  is  the  average  power  absorbed  per  unit  weight  of  the  biological  body. 
However,  there  appears  to  be  little  research  in  the  analysis  of  the  scat¬ 
tering  cross  section  of  arbitrarily-shaped  dielectric  and  biological  bodies. 

In  this  section  we  discuss  the  use  of  the  volume  integral  equation  to 
compute  the  scattering  cross  section  of  three-dimensional  arbitrarily-shaped 
dielectric  bodies  including  rectangular  and  I-shaped  boxes,  spheres,  finite 
circular  cy Linders,  and  simulated  birds. 

The  basic  volume  Integral  equation  has  been  discussed  in  detail  in 
Reference  o.  The  dielectric  body  can  be  replaced  by  an  equivalent  volume 
current  J  such  that 


where  is  the  anguLar  frequency.  E  is  the  electric  field, 


and  c  are 


the  complex  permittivity  in  free  space  and  the  dielectric  body,  respectively. 
Tiie  volume  integral  equation  in  terms  of  the  unknown  J  is 

/G(r,rf)  J  ^rf)  dv  -  — r  J  (r)  =  E^r),  (4  7) 

,  i  o  l  c-e  o  J 


G  (£,  r’)  =  -jwg(I  +  v  v  )  exp  (-jk  j  r-r  T 

k  4:»  j  r-r’  ! 


E  (r)  =  incident  electric  field  intensity, 


f  =  Principe 
Jv  singular 


pal  volume  integral  excluding  the 


gular  point  at  jr-r! 


I  =  unit  dyad 


The  solution  of  Equation  (47)  can  be  carried  out  by  the  method  of 
moments.  The  dielectric  body,  generally  heterogeneous,  is  divided  into 
rectangular  box  cells  and  the  equivalent  current  is  expanded  into  a  series 
of  pulse  functions,  each  of  which  is  uniform  in  one  cell  and  vanishes 
outside  the  cell.  The  Dirac-delta  function,  defined  at  the  center  of 
each  volume  cell,  is  used  as  the  weighting  function.  By  taking  a  scalar 
product  on  both  sides  of  equation  (47)  with  a  weighting  function  and  inte¬ 
grating  over  V,  we  generate  a  system  of  linear  equations  which  is  then 
solved  numerically  on  a  computer.  The  scattering  cross  section  is  then 
computed  in  terms  of  the  equivalent  current  J_  by  numerical  integration. 

Numerical  computations  have  been  conducted  for  dielectric  and  bio¬ 
logical  bodies  of  various  shapes  including  cubes,  cylinders,  spheres, 
rectangular  and  I-shaped  boxes,  and  simulated  birds.  Good  agreements  have 
been  observed  for  the  field  distribution  inside  the  dielectric  body  in 
comparison  with  the  data  from  Michigan  State  University  [  6,36,37].  For 
scattering  calculations,  the  only  data  available  in  the  literature  were 
for  spheres  and  finite  circular  cylinders.  The  present  calculation  showed 
correctly  the  sharp  resonance  behavior  of  the  back-scatter  cross  section 
as  a  function  of  frequency.  But  the  frequencies  of  resonance  were  about 
20  percent  lower  than  those  based  on  the  Mie  series  computation.  This 
discrepancy  could  be  due  to  the  reduced  apparent  size  of  the  sphere  in  the 
simulation  using  rectangular  cells.  Agreement  with  the  finite  cylinder  is 
good.  These  results  are  presented  in  detail  as  follows. 

A.  Scatterers  of  Simple  Shapes 

Figure  16  shows  the  calculated  back-scatter  cross  section  for  a  finite 
dielectric  cylinder  in  comparison  with  the  data  from  Richmond  [39].  Figure  17 
shows  the  geometry  of  a  rectangular  box  of  saline  water  and  the  way  the 
volume  cells  are  divided  and  numbered.  The  calculated  field  distribution 
is  displayed  in  Tables  2a,  2b  and  2c  for  the  x-component,  z-component  and 
total  field  of  the  electric  field  intensity.  They  are  in  good  agreement 
with  Michigan  State  data  [36].  Figure  18  shows  the  calculated  back-scatter 
cross  section  of  this  rectangular  box  of  saline  water  in  comparison  with 
the  measured  data  obtained  at  the  Georgia  Tech  compact  range. 
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TABLE  II 


ELECTRIC  FIELD  DISTRIBUTION  IN  THE  36-CELL  RECTANGULAR  BLOCK 


E  -DISTRIBUTION  IN  THE  36-CELL 
RECTANGULAR  BLOCK 


CELL 

NO. 

PRESENT 

CALCULATION 

MICHIGAN 

STATE  1 

CALCULATED 

MEASURED 

1 

.0510 

.0518 

.0518 

2 

.0600 

.0573 

.0523 

3 

.0952 

.0976 

.1632 

4 

.0816 

.0862 

.0862 

5 

.0883 

.0827 

.0832 

6 

.1040 

.1090 

.2355 

7 

.0870 

.0935 

.0935 

8 

.1225 

.1180 

.1180 

9 

.1355 

.1410 

.2869 

10 

.0292 

.0339 

.0459 

11 

.1008 

.1050 

.0878 

12 

.1611 

.1570 

.2017 

13 

.0930 

.1020 

.1258 

14 

.  1956 

.2020 

.2020 

1  5 

.  2 546 

.2480 

.1256 

16 

.  1345 

.1460 

.1460 

17 

.2267 

i 

I  .2330 

.  2330 

18 


.  3294 


.3220 


.5410 


TABLE  II  (Continued) 


E  “"DISTRIBUTION  IN  THE  36-CELL 
2  RECTANGULAR  BLOCK 


CELL 

PRESENT 

MICHIGAN 

STATE 

NO. 

CALCULATION 

1 

CALCULATED 

MEASURED 

1 

.2127 

.2140 

.2167 

2 

.1719 

.1710 

.1647 

3 

.1326 

.1310 

— 

4 

.0982 

.0963 

.0963 

3 

.0955 

.0927 

.0927 

6 

.0730 

.0698 

— 

7 

.0930 

.0991 

.0925 

8 

.0723 

.0783 

.0820 

9 

.0325 

.0364 

— 

10 

.2207 

.  2230 

.2230 

11 

.1355 

.1360 

.1456 

12 

.0554 

.0544 

— 

13 

.0666 

.0626 

.0626 

14 

.0330 

.03^8 

.0344 

15 

.0342 

.0348 

— 

16 

.0964 

.1010 

.1010 

mu 

.0932 

.0963 

.0963 

H 

i 

1  .0543 

.0563 

— 
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Figure  18, 
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Comparison  between  calculated  and  measured  RCS  of 
square  boxes  of  saline  water  shown  in  Figure  18. 
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the  1  vh.  ped  :  *>:  of  saline  water.  The  agreements  between  the  present 
calculation  and  that  at  Michigan  State  University  [36]  are  good.  Figure 
20  shows  the  comparison  between  the  measured  and  calculated  hack-scat  ter 
cross  section  data  generated  at  Georgia  Tech.  The  disagreement  could  be 
partially  due  to  the  acrylic  box  as  indicated  in  the  figure.  Good  agree¬ 
ments  with  Michigan  State  data  were  observed  also  for  three  other  cases, 
i::-luui:  c  a  cube  and  two  rectangular  cylinders,  which  w7ill  not  be  pre- 

Co  ut  at  ions  were  also  made  for  the  dielectric  sphere.  Figure  21 
*  ■.  *  tne  back-scat  ter  cross  section  of  a  dielectric  sphere  with  a  re- 
.  l  a-  p-.-rr  i  t  tivity  of  2.592.  The  results  deteriorate  as  Ka  becon- 

tit  r  th.iu  Figures  22  and  23  show  the  computed  results  for  a  dielect 

wui  w.tli  a  complex  dielectric  constant  of  29. 43-jO.  158  using  128-cell 
ut.u  -o.-ceii,  respectively.  The  results  are  rather  disappointing  when 
v::.:ai  ed  with  data  generated  by  Burr  and  Lo  [AO],  as  shown  in  Figures  2c 
and  -  ir.e  resonance  frequencies  were  shifted  by  20  percent  and  the 
resonance  teaks  art  off  by  30  percent.  This  failure  in  predicting  re¬ 
sonance  phenomena  in  a  dielectric  sphere  by  the  volume  integral  equation 
approach  is  in  contrast  to  the  high  accuracy  achieved  for  the  calculation 
of  conducting  spheres  by  a  surface  integral  equation  approach  [41,42j. 

The  results  for  the  sphere  were  presented  in  the  1980  IEEE  AP-S 
Symposium,  Both  Professors  Chen  and  Nyquist  of  the  Michigan  State 
University  attended  the  meeting  and  commented  that  they  had  also 
observed  similar  difficulties  with  the  sphere.  They  have  further 
noticed  that  if  the  incident  field  exp(-jkr)  is  decomposed  into  a  sine 
and  cosine  terms,  the  solution  corresponding  to  one  term  is  well  behaved 
and  that  for  the  other  is  ill-behaved. 


TABLE  111 


ELECTRIC  FIELD  DISTRIBUTION  IN  THE  128-CELL  BLOCK 


Ex-DISTRIBUTION  IN  THE  128-CELL  BLOCK 


CELL 

NO. 

PRESENT 

CALCULATION 

MICHIGAN  STATE 

CALCULATED 

MEASURED 

1 

.0545 

.0573 

— 

2 

.0342 

.0353 

— 

3 

.0656 

.0676 

— 

4 

.0279 

.0294 

— 

3 

.0906 

.0879 

— 

6 

.0899 

.0887 

— 

7 

.0787 

.0807 

— 

8 

.0557 

1  .0558 

— 

.1028 

.1071 

— 

10 

.0165 

.0173 

— 

11 

.1103 

.1057 

— 

12 

.1027 

.1018 

— 

13 

.07  75 

.0769 

.  0770 

14 

.0681 

i 

.0681 

.0681 

13 

.09b  7 

.1017 

.1017 

16 

.0393 

i 

.0386 

.0518 
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Ey~[ 

)ISTR  IBUTION  1 

N  THE  128- CELL  BLOCK 

CELL 

PRESENT 

MICHIGAN  STATE 

NO. 

CALCULATION 

CALCULATED 

MEASURED 

1 

.0228 

.0239 

— 

2 

.0482 

.0501 

— 

3  , 

.0467 

.0424 

— 

4 

.0698 

.0623 

— 

5  | 
| 

.0743 

.0681 

— 

1 

6 

.0497 

.0466 

— 

7 

.0309 

.0303 

— 

8 

.0867 

.0858 

— 

9 

.1254 

1  .1256 

10 

.1215 

.1234 

— 

11 

.0892 

1 

.0912 

— 

12 

.0490 

.0500 

i 

13 

.0197 

.0192 

.0243 

1 4 

.0450 

.0420 

.0316 

15 

.0496 

.0423 

.0423 

16  .0565  .0486  .0486 
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Ey-DISTR!BUT!ON  IN  THE  128-CELL  BLOCK 


CELL 


MICHIGAN  STATE 


NO. 

calculation 

calculated 

MEASURED 

17 

.0713 

.0671 

.0547 

18 

.0567 

.0544 

— 

19 

.0170 

.0162 

— 

20 

.0512 

.0485 

— 

21 

.0957 

.0937 

— 

22 

.1370 

.1382 

— 

23 

.1482 

.1506 

— 

24 

.1131 

.1136 

— 

25 

.0299 

.0296 

— 

26 

.0658 

.0641 

— 

27 

.0580 

.0544 

— 

28 

.0304 

.0262 

j 

— 

29 

.0164 

.0162 

.0194 

30 

.0375 

.0369 

.0371 

31 

.0355 

.0344 

.0344 

32 

.0165 

.0157 
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E  “DISTRIBUTION  IN  THE  128- CELL  BLOCK 


CELL 

NO. 

PRESENT 

CALCULATION 

|  MICHIGAN  STATE 

CALCULATED 

MEASURED 

1 

.05906 

.06208 

— 

2 

.05910 

.06123 

— 

3 

.08048 

.07979 

— 

4 

.07517 

.06709 

— 

5 

.11720 

.11119 

— 

6 

.10280 

. 10019 

— 

i 

.08449 

.(‘8627 

— 

8 

.10300 

.10235 

— 

9 

.16220 

.16500 

— 

10 

.12270 

.12460 

— 

11 

.14180 

.13961 

— 

12 

.11380 

.11342 

— 

13 

.  0  7  9  °  5 

.07926 

— 

14 

.08163 

.08001 

— 

15 

. 10870 

.11015 

— 

16 

.06883 

.06206 

39 
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Figure  20.  Comparison  between  calculated  and  measured  RCS  of  I-shaped 
boxes  of  saline  water. 
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B.  Scattering  Computation  for  1-foot  birds 


The  Green  Wingtail,  which  is  an  important  migrant  bird,  was  selected 
for  extensive  measurement  and  computation.  This  bird  is  typically  14.75 
to  15.50  inches  and  weighs  about  a  pound.  At  1  GHz  this  bird  is  about  1 
wavelength  in  length.  Figures  13  and  14  show  plaster-of-paris  models  for 
a  Green  Wingtail  in  sitting  and  flight  positions,  respectively. 

Four  simulation  models,  SB1  through  SB4,  were  generated  in  this  project 
For  the  first  three  models,  a  number  of  problems  developed  in  the  ex¬ 
perimental  work.  The  permittivity  was  too  high  or  too  low  or  not  uniform. 
These  difficulties  and  the  unsatisfactory  data  for  the  sphere  directed 
the  computation  toward  a  more  conservative  approach.  It  was  then  decided 
that  measurement  for  the  bird  should  be  started  yith  a  smaller  bird  exactly 
as  the  numerical  model  made  of  a  group  of  cubic  volume  cells.  With  this 
principle  in  mind,  SB4  was  fabricated.  Figure  24  shows  the  print-out  of 
the  cell  centers  for  the  side,  front  and  top  views  for  SB4.  Figure  26  shows 
the  geometry  of  the  coordinate  system  and  a  plane  wave  incident  in  the  x-z 
plane,  which  is  the  plane  of  symmetry  for  the  bird.  Figures  27  and  28  show 
fair  agreement  between  the  computed  and  measured  results.  It  is  noticed 
that  the  agreement  is  good  as  long  as  there  is  no  sharp  variation  in  the 
pattern.  This  difficulty  in  predicting  sharp  field  variation  and  resonance 
phenomena  may  explain  the  use  of  specific  absorption  rate  (SAR)  in  dealing 
with  dielectric  scattering  problems  [38]. 
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BACK  SCATTERING  CROSS  -  SECTION  (<1B) 


Figure  do.  Hack,  scattering  crosa  section  oi  bird  SH4  versus  elevation 
angle  ■■  with  vertical  polarisation  as  shown  in  Figures  .  ’ 


SECTION  V 

IMPROVEMENTS  OF  COMPUTER  ALGORITHM 


A  number  of  modifications  have  been  made  on  the  existing  Georgia 
Tech  volume  integral  equation  algorithm  to  make  it  more  efficient  for 
the  computation  of  dielectric  scatterers,  Major  improvements  included: 

(1)  the  reduction  of  execution  time  and  central  memory  requirement  by 
50%  by  using  symmetrical  matrices,  and  (2)  the  reduction  of  execution 
time  and  central  memory  requirement  by  75%  for  scattering  problems  with 
one-plane  symmetry  and  by  87.5%  for  scatterers  with  two-plane  symmetry. 

In  addition,  the  banded  matrix  [43,44]  and  virtual  memory  [45]  techniques 
have  been  implemented  successfully  for  small  scatterers  and  a  new  local- 
file  manipulation  technique  has  been  explored. 

A.  Symmetrical  Matrices 

It  can  be  shown  that  the  matrix  of  the  volume  integral  equation  al¬ 
gorithm  is  symmetrical  if  the  following  conditions  are  satisfied. 

(1)  the  dielectric  body  is  homogeneous  and  has  a 
constant  permittivity  » 

(2)  the  volume  cells  are  equal  in  size  , 

(3)  the  volume  cells  are  identical  in  shape. 

The  third  condition  is  not  critical  and  can  often  be  ignored.  To  prove 
this,  one  can  examine  the  following  matrix  elements  in  the  algorithm 
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where 


3ju:'U(r,)-Le]  °n 


k,n  =  1,2,3  or  x,y ,  z 


i,p  =  1,2,. ...L;  the  index  number  of  volume  cell, 
xk  ®  the  Kronecker  delta  function. 


®  a  pulse  function  being  unity  in  the  2th  volume 
cell  associated  with  k  (which  designates  x,y  or 
z  and  is  merely  a  dummy  index  number  of  no 
consequence  to  the  integration). 


Since  the  pulse  function  B  restricts  the  domain  of  integration  to  the 

p 

unit  volume,  cell  1  in  Equation  (4S)  and  cell  p  in  Equation  (49), 


r  —  r  and 


are  equal.  Thus  /  ,  and  Z*  are  eaua!  under  the  tnree 
*  p  £  n 


t  c  r  n  s v rime tneai 


itrix,  rhere  are  standard  subroutines  to  handle 


solution.  The  improved  Georgia  Tech  algorithm  is  given  an 


'  B P WSM"  to  denot  ( 


lioiogical-Pla  n  o  -  V.Ta  ve  -  S v  mm  e  t.  r  1  c  a  1  Matrix". 


i  n  i  a  j  g  o  r  1 1  hm  has  b  o  o n  s  u c cos  •"  u 


n ur.he  r  i,-  J  o n  e  s  v:  i  t 1 


.l a  e . j :n p liter  . .  nr  rai  mo m:  r v  a n u 


:  u  1 1  o  n  1 1  mi' . 


Two  improved  o'  gori  thms ,  B!iPV2  and  BHPW4 ,  have  been  completed  tc 


r  e  s  u  i.  l  i.  n  g  r  - 


.it  taring  r rob i e rr,  with  one-oiane  and  two- plane  s*.*mmet ries  with 
redu  t .  or  /Oa.  ana  b/.oi  in  computer  time  ana  central  memory 


./e-tiveiv.  These  are  discussed  separately  as  i widows. 

i)  scat ter iin  problems  with  one-pi ane  symmetry  ibtiPWlj 

I'nitn  the  direction  oi  propagation  ot  an  incident  p.ane  wave 
in  tne  plane  of  symmetry  oi  a  conducting  scatterer,  it  is  re- 
,  ..  ;cud  Ltiat  some  symmetric  behavior  must  exist  in  the  induces 

ou  the  surface  ut  the  scatterer.  Without  loss  oi  generality, 
- 1  ■  i : i  coordinates  can  be  set  up  so  that  the  plane  ot  symmetry 

*ith  the  XZ  plane  as  shown  in  Figure  db.  The  no larization 
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urt*  2  6.  A  symmetrical  si  at  r  i  i  I  um  i  11.1 1  • i>i  *  .  i  cl  in.-  wti,f  in;i<hnt  in  the  plane  oi  symmetry 


of  the  incident  wave  is  assumed  to  be  either  parallel  or 
perpendicular  to  the  XZ  plane.  Arbitrary  polarization  can 
be  decomposed  into  two  components,  one  parallel  and  the  other 
perpendicular  to  the  XZ  plane.  The  overall  scattering  problem 
can  then  be  treated  by  superposing  the  fields  due  to  these  two 
component  incident  fields. 

At  two  symmetrical  points  £  and  £  +  L/2  in  Figure  28,  the 
components  of  the  induced  currents  exhibit  the  following  relation¬ 
ships 


Jx  =  Tx 

£  J£  +  L/2  ’ 


Ty  =-Ty 


£  +  L/2  , 


jz  =  jz 

£  £  +  L/2  , 


(50) 


when  the  incident  E  is  parallel  to  the  z-axis.  In  Equation  (50 ) > 
J*  denotes,  for  example,  the  x  component  of  the  induced  current  in 
the  £th  cell. 

When  the  incident  E1  is  parallel  to  the  y-axis,  the  induced 
current  on  the  scatterer  has  the  following  property 


JX  =  - JX 

J£  J i  +  L/2  , 


jy,  =  jy 


l  +  L/2  ’ 


JZ  =“JZ 

£  J£  +  L/2  . 


(51) 


The  matrix  equation  to  be  solved  is 


1  3  1  0 

y  v  jk  z*?  =  vn 

v*=i  k=i  ^  pk  p 


n=l ,2,3 


p=l, ...L. 


Since  the  excitation  is  symmetrical,  we  have 


vn  =  vn 

l  \  +  L/2 


Substitution  of  Equations  (50)  and  (53)  into  Equation  (52) 

yields 


E  Jk  !  Z"“  +  Z(L/2‘K>n  k+1]  =  v” 

k=l  c  P 


n=l,2, J 

p=l,2 - L/2 

for  ItE  =  z 


Substitution  of  Equation  (51)  and  (53)  into  Equation  (52) 
yields 


L/2  3 


Jt.^n  ^  _ (L/ 2+i)n 


X  E  JI[zpk  +  zpk 


n=i  k-=i 


n=l ,2 , 3 
p=l , 2 , L/2 
for  =  y 


•  nations  (54) 


•  L  '•-*  ~  p  1  ci  LI  L"  S 


id  (  ■•’»  ^how  that  the  number 

.  >  '  ♦  . .  -  '  n 


r  >n-  has 


Figure  29  shows  a  scattering  problem  symmetrical  with 

respect  to  the  x-z  and  y-z  planes.  The  directive  of  propagation 

of  a  plane  wave  is  assumed  to  be  parallel  to  the  z-axis .  l%hen 
a 

L  -  '■*,  we  nave 
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in  addition,  the  excitation  at  cells  1,2,1  and  4  is  identical. 
The  number  of  Equations  in  Equation  (34)  can  therefore  be  reduced 
to  one-quarter  of  its  original  size  as  follows 


L/4  3  ,  |  ,  ,.il  f(L/4+t)«  (L/2+«)n 

E  E  J«  |zpt  +  <-»  IV  +v 

7=1  k=l 


(3/A  L+H)n 


n  =  1,2,3 
P  =  1,2,  . 


.,  L/4 


when  E  =  z  . 


Substitution  of  Equations  (57)  into  Equation  (52)  yields 


L  r  Ji  (-Dk(1“5lk)  z(k/4+£)n 

i-j i  k^i  M  PK  pk 


+  (-i)k(1-6Zk>  z  <L/2+2)n 

pk 


+  (_l)k(l-63k)  z  L4'£)n  ?  _  vn 


n  -  1,2,3 
P  =  1*2,  . 


.  . ,  L/4 


E  =  y.  In  Equations  (59)  is  the  Kronecker  delta. 


C.  Banded  Matrix  Techniques 

The  banded  matrix  technique  has  been  previously  employed  by  Ferguson, 
et  al.  [43]  and  Balestri,  et  al.  [44]  in  the  scattering  and  radiation  of 
thin-wire  structures.  These  authors  have  demonstrated  that  the  banded 
matrix  technique  can  reduce  the  computer  execution  time  in  the  computation 
of  thin-wire  scattering  and  radiation  involving  small  matrices.  They  also 
demonstrated  that  wire-grid  problems  involving  more  wire  segments  than  can 
be  managed  in  the  computer  central  memory  can  be  solved  by  the  banded  matrix 
technique  after  proper  numbering  of  wire  segments.  The  segments  are 


77 


numbered  such  that  the  difference  between  segment  numbers  for  all 
neighboring  segment  pairs  is  small  as  compared  with  the  total  number 
of  segments.  In  this  case  the  large  matrix  elements  are  kept  close 
to  the  principal  diagonal  of  the  matrix. 

The  basic  banded  matrix  technique  has  been  applied  to  the  Volume 
Integral  Equation  algorithm.  However,  there  are  three  unknowns  in  the 
nth  cell  in  the  Volume  Integral  Equation  approach,  while  there  is  only 
one  unknown  in  the  nth  segment  in  the  wire  algorithm.  Fortunately  the 
matrix  generated  in  the  Volume  Integral  Equation  algorithm  has  a  tendency 
to  be  banded.  Figure  31  shows  the  matrix  for  the  case  of  a  simple  cylinder 
formed  with  a  linear  array  of  volume  cells.  It  can  be  seen  that  a  diagonal 
band  extended  to  one  third  of  the  columns  and  the  rows  must  be  included  in 
order  to  include  all  the  nonzero  elements.  It  is  also  noted  that  outside 
this  band  all  the  matrix  elements  are  zero.  This  phenomenon  is  due  to  the 
lack  of  coupling  between  orthogonal  components  of  the  electromagnetic 
source  and  field.  For  an  object  of  more  complex  geometry,  such  as  a  sphere, 
elements  throughout  the  entire  matrix  can  be  nonzero  except  for  those  re¬ 
lated  to  the  coupling  between  orthogonal  components  in  the  self  cell. 

Figure  32  shows  the  matrix  for  the  case  of  a  prolate  spheroid  of  12  cells, 
which  is  weakly  banded  with  some  nonzero  elements  away  from  the  diagonal 
band . 

We  now  define  the  normalized  width  of  the  diagonal  band  as 

Normalized  Width  of  diagonal  band  =  [  (number  of  rows  in  band) 

+  (number  of  columns  in  band)]  [(total  number  of  rows)  x  2]  (60) 

Numerical  tests  have  been  conducted  to  explore  the  convergence  of  the 
solutions  as  a  function  of  the  normalized  width  of  the  diagonal  band. 

Figure  33  shows  that  the  error  in  the  solution  is  reduced  as  the  width  of 
the  diagonal  band  is  increased.  The  error  falls  to  near  zero  when  the 
normalized  width  of  the  diagonal  band  is  only  0.6.  Note  that  for  a  general 
matrix  the  error  reaches  zero  when  the  normalized  width  of  the  diagonal 
band  approaches  unity.  It  is  also  noteworthy  that  the  error  is  only  5% 
when  the  normalized  width  of  the  diagonal  is  0.13  or  larger.  Thus,  a 
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Note:  0  The  size  of  the  shaded  area  in 
the  square  cell  is  proportional 
to  the  magnitude  of  the  matrix 
element. 


Figure  i(j .  A  display  of  the  magnitude  of  matrix  elements  in  a  VIE 

matrix  for  an  8-cell  dielectric  cylinder  which  is  strongly 
banded . 


y  n 


D.  Virtual  Memory  Technique 

The  virtual  memory  is  a  software  technique,  in  contrast  to  the 
extended-memory  hardware  capability,  to  store  data  in  a  digital  computer 
for  rapid  and  efficient  access  in  computer-aided  numerical  analyses  to 
overcome  the  limitations  imposed  by  the  size  of  the  computer  central 
memory.  The  basic  algorithm  had  been  developed  by  Carbrev  [45]  for  real¬ 
valued  data,  and  was  employed  in  the  present  research  to  expand  the 
capability  of  the  Volume  Integral  Equation  algorithm  to  handle  large 
dielectric  scatterers.  When  using  exclusively  the  computer  central 
memory  to  handle  the  matrix  in  the  computation,  the  CDC  CYBER-74  computer 
at  Georgia  Tech  can  only  deal  with  matrices  of  about  38,300  complex 

elements,  or  about  65  volume  cells.  The  use  of  the  virtual  memory  tech- 

5 

nique  can  potentially  make  it  possible  to  handle  a  matrix  with  4  x  10 
complex  elements,  or  660  volume  cells,  which  is  about  ten  times  the  size 
of  those  limited  by  the  computer  central  memory.  A  serious  disadvantage 
of  the  virtual  memory  technique  is  its  extremely  large  execution  time, 
often  ten  times  more  than  methods  using  central  core  memory  alone. 

The  tasks  involved  in  implementing  the  virtual  memory  technique  in 
the  Volume  Integral  Equation  algorithm  are  twofold;  the  Carbrey  algorithm 
must  be  extended  to  handle  complex  data  and  be  integrated  into  the  process 
of  the  moment  method  solution.  Both  of  these  difficulties  have  been  overcome 
and  successfully  tested  for  small  and  medium  scatterers.  For  large 
scatterers  considerably  exceeding  the  central  memory,  the  computational 
efficiency  is  low  and  needs  to  be  improved. 

There  are  two  steps  involved  in  the  numerical  solution  of  a  system 
of  linear  equations.  First  the  matrix  elements  involved  must  be  computed 
and  stored  for  easy  access.  Secondly,  the  matrix  equation  must  be  solved 
by  a  certain  process  using  either  the  central  memory  alone  or  the  virtual 
memory,  which  uses  both  the  central  memory  and  disk  memory.  Both  steps 
need  large  computer  execution  time  and  the  key  issue  is  to  reduce  the 
computer  time  to  a  level  acceptable  for  practical  computations. 
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F. _ Local-F lie  Jlan Ip u lat  i  o n  l  e eh a iy[ u e 

As  has  been  :u»ud  in  the  . .receding  section  ,  the  virtual  memory 
'.eehnique  is  c  r:  •  :  a  :  c  iarcr-  me  tcice.-  'cvopu  the  capacity 

tile  computer  ■_  ‘iiir..; .  at  n  •:  :  udiv.^. .. ,:ed  D  y  the  c:^a.  :i;dv 

ton  -  execution  time  ..eeJeti .  For  example,  15U,000  CPU  sec  is  recuired 
in  a  computer  run  involving  the  solution  of  a  48b  x  487  matrix.  Further¬ 
more,  tile  CPU  time  required  increases  astronomically  with  the  number  of 
matrix  elements  when  it  exceeds  the  capacity  of  the  computer  central 
memory.  In  order  to  overcome  this  difficulty,  a  highly  efficient 
’’local- f  i  1  l.  manipulat  ion”  technique  was  developed  at  Georgia  Tech, 
i his  technique,  which  iias  been  tested  only  for  snail  and  medium-sized 
matrices,  utilizes  both  the  discs  and  the  central  memory  for  the  storage 

similar  lo  the  virtual  memory  technique,  the  local-file  manipulation 
net nod  is  a  software  technique  and  it  has  been  written  at  Georgia  Tech  in 
Fortran  .  and  COMPASS  languages:  the  latter  being  a  machine  language  used 
ea  CjC  CvUeir  7  u  ihev  ke*-  elements  of  the  program  include  twent v -seven 
s ub rout  i ne s  ,  e i eve n  o f  wii iuii  are  write en  in  COMPAF S .  Tlie.se  subroutines 
were  imp  1 counted  in  the  dielectric  scattering  calculation  urogram  used  in 
this  research  project-  ihev  can.  be  divided  into  three  groups  according 
to  ihoir  functions  as  illustrated  in  Figure  3.3  and  Table  IV.  Appendix  I 
list.-  all  these  subroutines  and  describe  their  functions. 

The  kev  features  of  the  local- file  manipulation  technique  are  in 
tiie  efficient  wavs  of  storing  and  recalling  the  matrix  elements  as  well 
as  the  specific  method  of  solving  the  system  of  linear  equations  by  the 
1 owe r -u p pe r  de compo s i t i on  method . 

All  the  matrix  elements  are  stored  in  the  present  program  in  four 
local  files,  each  of  which  is  allocated  a  certain  disc  space.  The  matrix 
elements  are  first  divided  into,  say,  three  portions  as  shown  in  Figure  34. 
Elements  in  each  of  the  corresponding  rows  of  Lhe  three  partitions  are  computed 
and  stored  in  small  arrays  in  the  core  memory.  As  soon  as  the  three  rows 
have  been  fuLiv  generated,  these  three  rows  are  transferred  from  core 
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memory  to  local  files  in  binary  form  with  a  "WRITE"  statement.  The 
process  is  repeated  for  N  times  until  the  complete  matrix  has  been 
generated  and  stored  in  the  three  local  files.  This  storage  process 
is  highly  efficient  in  comparison  with  the  conventional  virtual  memory 
technique  because  matrix  elements  of  entire  rows  are  read  simultaneously 
into  local  files. 

The  solution  of  the  matrix  or  the  system  of  linear  equations  in  the 
present  local-file  manipulation  technique  uses  the  lower-upper  decomposition 
method  described  in  Appendix  II.  There  are  two  important  features  in 
this  technique;  the  restorage  of  the  matrix  elements  based  on  their 
frequencies  of  impending  usage  and  the  solution  of  the  matrix  by  the 
lower-upper  decomposition  method. 

In  the  restorage  of  the  matrix  elements  a  priority  table  is  generated 
and  stored  in  core  memory  to  assign  the  "priority"  level  for  each.  Three 
priority  levels,  high,  last  and  low,  are  assigned  to  the  rows.  The  "high" 
priority  is  assigned  to  the  row  for  which  immediate  and  frequent  access 
is  needed.  The  local-file  manipulation  algorithm  has  been  tested  on  the 
Georgia  Tech  CDC  Cyber  74  computer  with  good  results  for  small  and  medium 
sized  matrices.  Figure  35  shows  a  comparison  for  the  time  required  for  the 
generation  and  storage  of  the  matrix  elements  between  the  virtual  memory 
algorithm  of  Carbrey  and  the  present  local-file  manipulation  algorithm. 

This  greatly  improved  efficiency  in  the  storage  process  is  probably  due 
to  the  row  storage  method  discussed  previously.  The  matrix  solution 
time  required  in  the  local-file  manipulation  technique  is  also  much  less 
time  consuming  than  the  Carbrey  algorithm  as  shown  in  Figure  36.  These 
comparisons  are  made  only  for  small  matrices  and  it  remains  to  be  seen 
whether  similar  comparisons  hold  for  very  large  matrices.  Unfortunately 
further  study  for  the  case  of  large  matrices  cannot  be  conducted  within 
the  time  and  resources  of  the  present  research  program.  It  appears,  how¬ 
ever,  that  the  local-file  manipulation  technique  will  prove  to  be  con¬ 
siderably  superior  to  the  Carbrey  virtual  memory  technique  when  applied 
to  large  matrices.  A  major  advantage  of  the  lower-upper  decomposition 
process  is  that  one  does  not  have  to  repeat  the  matrix  element  generation 
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i  or  dif  ferent  -■:<  ji  tuitions  and  ther^  r-...-ro  save.;  large  C.pr  when 

tire  angle  and  polarization  of  illumi"  nioi:  *'•  changed.  A  d  «  c.^v-j  ntagt- 
of  this  method  is  "hat  it  requ  L  r*-.s  .vie,-  as  much  -tora-e  spa.;-  as 
other  methods  sac1'  as  the  (hiu  ssian-  - 1  Lmi  ea  tion  method  because  of  the 
need  to  store  both  the  upper  and  lower  matrices  [L]  and  [ U ]  .  However, 
this  disadvantage  is  not  of  significant  consequence  when  it  is  used  in 
the  local-file  manipulation  method.  Since  the  matrices  are  stored  in  the 
virtually  unlimited  disc  space,  the  efficiency  of  storage  and  access  to 
tiie  matrices  are  more  important  than  the  size  of  the  storage  space 
required.  In  addition,  the  enlarged  stord^.;  requirement  in  this  method 
can  be  reduced  by  the  band-storage-mode  technique  [ 4 6 J . 
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SECTION  VI 


METHODOLOGY  FOR  THE  MODELING  OF  THE  SCATTERING 
! '  A  i  LOCK  OF  FLYING  Ti.vDS 

la  radar  applications,  the  electromagnetic  scattering  from  a  flock 
of  flying  birds  interferes  with  the  detection  and  identification  of 
objects,  but  it  can  also  be  used  to  track  bird  migratory  patterns  [47,48]. 
Radar  return  from  birds  is  often  observed  as  a  large  dot  angel,  which 
occurs  at  all  microwave  frequencies  on  virtually  every  type  of  radar.  One 
of  the  most  outstanding  features  of  the  scattering  from  a  flock  of  birds, 
and  hence  the  radar  angel,  is  its  irregularity;  with  variations  of  as 
much  as  three  orders  of  magnitude  [49J.  This  irregularity  is  due  to  that 
of  the  individual  Dird  [2$]  and  the  density  and  formation  of  the  flock. 

Although  the  lack  of  repeatability  and  regularity  in  the  scattering 
measurements  for  birds  appears  to  fit  a  statistical  model,  several  aspects 
of  the  problem  are  highly  deterministic.  The  formation  of  a  specific 
flock  is  quite  orderly  and  therefore  wili  probably  be  more  accurately 
modeled  deterministically.  The  individual  birds  can  be  better  modeled 
statistically  because  of  the  variations  in  size,  shape,  and  movements 
among  them. 

Let  o  G j  ,  j )  be  a  random  variable  which  represents  the  back-scatter 
cross-sec t ion  irom  a  specific  ivpe  of  bird  as  shown  in  Figure  37,  the 
total  radar  cross-section  o  can  be  approximated  by 
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where  r  is  a  unit  vector  parallel  to  the  line  of  sight  between  the  radar 

and  a  reference  bird  in  the  flock,  r  is  the  distance  from  the  reference  bird 

-n 

to  a  bird  designated  as  number  n,  N  is  the  total  number  of  birds  in  the 
flock.  Let  the  superscript  denote  ’’average",  we  have  the  expected  total 
KGS  of  the  flock  as 
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(62) 


v  =  .*)  V 


where  it  is  assumed  that  the  formation  of  the  flock  is  determines tic , 
at  least  for  the  time  interval  under  consideration. 

The  time  variation  of  the  radar  return  from  a  flock  of  flying  birds 
is  primarily  due  to  the  wing  flapping  and  the  flight  path  movement.  The 
effect  of  this  variation  appears  as  radar  scintillation  and  Doppler 
frequency  shifts.  Modeling  of  these  phenomena  can  be  based  on  the 
statistical  characteristics  of  the  random  variable  a  ( t , y )  pertaining  to  a 
single  bird. 

Depending  on  the  available  data  on  the  single  bird,  there  are  several 

ways  to  model  the  scattering  characteristics  of  a  flock  of  birds  by  using 

equations  foi)  and  (62).  Equation  f61)  can  be  used  to  generate  a  Monte- 

Carlo  simulation  with  the  statistical  parameters  of  the  single  bird.  If 

it  is  desired  to  include  the  effect  of  the  arrav  formation,  r  can  also 

~n 

be  considered  as  a  random  variable  and  do  included  in  the  model..  If  we 
are  only  interested  in  the  time  average  scattering  cross-section.  Equation 
(b2)  is  convenient  to  use  as  long  as  the  average  scattering  cross-section 
oi.  a  single  bird  is  known.  The  singLe  bird  data  can  be  obtained  by 
computation  and  measurement  techniques  discussed  in  this  report.  The 
techniques  oi  modeling  a  flock  ol  birds  from  scattering  data  of  a  single 
bird  is  similar  in  many  ways  to  that  of  chaff  cloud  scattering  problem, 
which  has  been  recently  investigated  at  Georgia  lech  [50-32].  No  numerical 
modeling  was  carried  out  in  the  research  program  because  of  the  lack  of 
quality  data  which  can  be  used  to  compare  with  the  model. 


SECTION  VII 

CONCLUSIONS  AND  RECOMMENDATIONS 


Research  has  been  conducted  in  the  analysis  and  measurements  of 
three-dimensional  arbitrarily-shaped  heterogeneous  dielectric  and 
biological  bodies.  The  discrepancies  in  the  literature  regarding  the 
singularity  of  the  electric  dyadic  Green* s  functions  were  resolved. 

The  discrepancies  were  centered  at  the  singularity  of  the  rectangular 
cavity.  It  was  shown  in  this  report  that  the  apparent  discrepancies  do 
not  exist  and  a  unified  and  consistent  view  was  presented. 

Compact  range  scattering  measurements  were  successfully  conducted 
at  1  GHz.  Techniques  in  fabricating  simulation  models  using  the  "Super- 
stuff"  were  investigated  and  several  1-foot  birds  were  made.  Extensive 
numerical  analysis  was  carried  out  for  dielectric  scatterers  of  various 
shapes  including  cylinders,  rectangular  blocks,  I-shaped  blocks,  spheres, 
and  a  1-foot  bird.  The  accuracy  of  these  computations  was  good  except  for 
the  resonant  sphere,  for  which  the  resonance  frequencies  were  shifted  by 
about  20  percent. 

Various  numerical  techniques  have  been  investigated.  Computer  central 
memory  requirements  and  execution  time  were  reduced  by  50  percent  with  the 
symmetrical  matrix  technique  and  by  75  percent  for  scatterers  of  one-plane 
symmetry.  Banded  matrix  and  virtual  memory  techniques  have  been  implemented 
in  the  Volume  Integral  Equation  algorithm  and  tested  successfully  for  small 
scatterers.  A  new  local-file  manipulation  technique  for  handling  matrices 
larger  than  core  capacity  was  explored  and  the  results  for  small  matrices 
showed  that  it  is  potentially  much  more  efficient  than  the  virtual  memorv 
algorithm  developed  by  Carbrev. 

It  is  recommended  that  the  computational  techniques  developed  in  this 
research  program  be  extended  to  the  analysis  of  larger  lossy  dielectric 
bodies.  In  addition,  the  Surface  Integral  Equation  approach  should  be 
expLored  to  see  whether  its  numerical  convergence  is  more  rapid  than  the 
Volume  Integral  Equation  approach.  The  derivation  of  an  exact  soLution 


q-t 


for  the  dielectric  prolate  spheroid  appears  to  be  a  feasible  research 
subject  which  should  lead  to  accurate  and  useful  data  for  dielectric 
scattering  problems. 
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Sixteen  of  the  subroutines  are  written  in  FORTRAN  V  and  eleven  of 
them  are  in  the  COMPASS  machine  language. 

A.  Subroutines  Written  in  FORTRAN  V 

1)  MMCDCOM  -  Decomposes  a  complex  matrix  to  lower  and  upper 
submatrices . 

2)  MMCFBS  -  Solves  the  decomposed  matrix  using  forward  and 
backward  substitution. 

3)  MMCRNRM  -  Normalizes  a  row  in  a  matrix  by  dividing  each 
(except  the  first)  element  of  the  row  by  its  first  element. 

The  first  element  of  the  new  row  is  then  replaced  by  its 
reciprocal . 

4)  MMCRSUB  -  Performs  complex  row-subtraction  and  the  pivoting 
operation  in  the  lower-upper  decomposition  method. 

5)  MM1NIT  -  Pre-sets  matrix  storage  allocation  among  several 
local-files,  allocates  local-file  addresses  and  allocates 
buffers  for  the  first  few  rows.  This  is  an  initialization 
of  parameters,  which  may  be  altered  during  the  execution 
of  the  program. 

b)  MMSTAT  -  Provides  statistical  information  regarding  the  status 
and  the  efficiency  of  the  core  and  buffer  storage  allocation 
as  well  as  other  program  execution  characteristics. 

7)  MMGBUF  -  Searches  for  and  provides  a  free  buffer  upon  request. 

8)  MMINBND  -  Assigns  the  priority  cor  a  row  based  on  whether 
it  is  of  current  interest  or  not. 

9)  MMFSET  -  Sets  up  a  file-environment  table  for  an  I/O  operation 
and  starts  the  I/O  operation. 

10)  MMLOCK  -  Assigns  buffers  and  insures  that  the  buffer  is  available 
for  usage. 

11)  MMGTRA  -  Gets  the  disc  address  of  a  matrix-row  from  the  address 
tables  such  as  the  row  map,  buffer  maps,  FET,  etc. 

12)  MMSINBD  -  Sets  a  map  of  in-bound/out-of-bound  rows  for  later 
use  by  MMINBND. 

13)  MMGRDW  -  Gets  a  matrix  row  from  disc  and  stores  it  into  the  core 
through  the  buffer. 


* 


14)  MMSROW  -  Gets  a  matrix-row  from  the  core  and  stores  it 
into  the  disc  through  the  buffer.  This  is  the  reverse 
process  of  MMGROW. 

15)  MMSLIM  -  Determines  and  updates  the  priority  table  which 
assigns  the  matrix  rows  to  be  stored  into  the  core  memory. 

16)  MMBAHED  -  Gets  matrix  rows  to  be  used  soon  from  the  disc 
and  stores  them  into  the  buffer.  This  is  a  portion  of 
MMGROW  and  is  not  needed  for  MMSROW. 

B.  Subroutines  Written  in  COMPASS 

1)  MMALLOC  -  Allocates  central  memory  space  for  the  buffers 

to  be  used;  Sets  up  parameters  such  as  the  number  of  local - 
files  to  hold  the  generated  matrix,  the  number  of  rows  per 
super-row,  and  the  number  of  simultaneous  I/O  requests  that 
are  allowed. 

2)  MMLINK  -  Links  subroutines  which  allocate  buffers,  set  file 
environment  tables,  row-maps,  buffer-maps,  and  subroutines 
indirectly  calling  other  subroutines. 

3)  MMSUPRT  -  Contains  a  collection  of  subroutines  to  direct 
the  computer  to  perform  several  simultaneous  operations  to 
take  advantage  of  the  overlapping  mode  of  execution  of  the 
Cyber- 74 . 

Y)  MMCIO  -  Invokes  MCIOM  to  process  FET  (file  environment  table). 

FKT  is  a  file-table  containing  information  such  as  the  address 
of  the  randomly-accessed  records  of  a  particular  file.  "CIO" 
stands  tor  confined  L/0  operations  and  is  used  to  perform  the 
I/O  operation  in  COMPASS.  In  FORTRAN ,  COBOL  and  other  high 
level  languages  the  I/O  is  performed  automat icaliy . 

))  MMMSG  -  DispLays  messages  on  the  computer-operator  console 
and  opt  ion  illy  display.,  messages  in  Lhe  daylile. 

o)  (  ol)  -  converts  binary  displav  to  decimal  display. 

/)  MMKCL  -  M- *n :  t  ors  the  1/0  process  to  insure  tiiuL  Lhe  i/0 
I  pe  n;  it'll  1  enmp  1  eled  . 

o)  MMCOPY  -  Copies  data  iron  one  array  to  the  other. 

u)  MMKBl'sY  -  Checks  i)  an  1,0  t  or  a  particular  super-row  is 

Ij  )  MMSPLIM  -  Sets  the  PPL'  limits.  fA  periplieral  processor  unit 
i P Pi J >  is  a  small  processor  whicn  provides  communication  paths 
between  the  central  processor  and  the  individual  peripheral  devices 
such  as  the  disc  unit.  in  this  subroutine  the  user  can  .Limit 
or  expand  the  number  of  simultaneously  executed  I/O  requests 
depending  on  the  computer’s  load,  time  of  day,  etc. 

11)  MMUPACK  -  Unpacks  the  row-map  and  buffer-map.  Keeps  track 

on  a  table  of  l lie-numbers  corresponding  to  row,  super  rows, 
and  buffers.  Checks  ii  a  butler  cannot  be  released  lor  a 
part  icul  .ii  row  . 
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MATRIX  SOLUTION  BY  THE  LOWER-UPPER  DECOMPOSITION  METHOD 


The  lower-upper  decomposition  method  for  the  solution  of  a  matrix 
equation  is  based  on  the  following  theorem  , 

Theorem.  Let  [A]  denote  a  NxN  square  matrix  whose  elements  are  denoted 
by  a  ,  where  i  and  i  refer  to  the  row  and  column  of  the  element.  Let 

ij 

[A^J  denote  a  kxk  square  matrix  which  is  a  left-upper  submatrix  of  [A] 
defined  as 


[AjJ  =  [a..];  i , j  =  1,  ...k.  (AI-1) 

If  all  the  submatrices  [A^J,  k=l,...N,  are  nonsingular,  then  f A ]  can  be 
decomposed  uniquely  into  the  following  form 

IA]  =  [L]  [UJ  (AII-2) 

where  all  the  three  matrices  are  all  of  dimensions  NxN,  and  [ L ]  and 
L U J  are  referred  to  as  the  lower  and  upper  triangular  matrices  of  the 
following  form 


L= 


1  0  0, 

i21  1  0, 

*nl  ln2  ' 


L11  U12 .  Uln 

J  .  U22 .  ^2n 

U  0  . U  . 

nn 


(AII-3) 


(AII-4) 
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‘cnnsferring  [A]  into  :;ie  product  of  twe  such  matric un,  one  can 


sclve  the  matrix  equation  b>  f  o rward  and  backward 


. :  ru  L  LOi« 


in*,  d  iat 

er. 

W C  [,v'i  :  i  S .  1 J W  O'.' 

•.  [L;  I'ld 

— 1 

P 

H ' 

U1  2 

u13 

[A]  = 

*21ull 

fc21U12+u22 

Jt2iui3+u; 

SlUll 

^  3 1 U1 2^^  32  u2  2 

l3iun+i: 

[AII-5] 


11 

a12 

a13 

21 

a22 

a23 

nL 

a  „ 
n2 

an3 

By  equating  the  individual  elements  in  the  equation  above,  we 
obtai n : 


a)  For  the  first  row  of  U 


u  =a 

ij  ij 


b)  For  the  first  column  of  L 


j-1,2. 


iil  Uil/u. 


1=1,2, 


c)  For  the  second  row  of  U 


u  ,  ,=a  .  .  -1 , .  u.  . 
2j  2j  21  lj 


j = 2 ,  J , 


d)  For  the  second  column  of  L 

1i2_  (uia~iilu12^/u22 


i=2 , 3 , 


''he  third  mw  -j  f 


All.  the  matrix  elements  t,  .  and  u .  .  can  then  be  obtained  in 

aj  J-J 

sequence  according  to  the  set  of  Liquations  (All-6)  .  Note  that  the 

simple  process  of  determining  the  unknown  e..  and  u.;  involves  only 

st ra i gh t forward  substitution  of  the  known  a. .  and  tne  previously 

calculated  «.  .  .  and  u.  .  following  the  order  specified  in  Lquation 
:  )  hj 

Note  that  hqua l  ion  ( A 1  i. -  6 )  can  be  denoted  in  general  as 


,  u,_  j*P,pfl,...N 

d.i  /  *vj 


=  (a. 

1  q  1  q 


E 

k«l 


-  U  )/u  i=q+l ,  .  .  .  N  (  at  j  —  y  ) 

ix  kq  qq  M  votL  / ) 


Let  us  now  use  this  lower-upper  decomposition  method  to  solve  the 
i  o  iviwing  matrix  equation 


l A J  fXj  =  fh  | 

Appiving  liquation  (AI-2)  to  (Al-8) ,  we  have 

h  i  ;rj  (x i  =  [b] 


(AII-8) 


(All-9) 


;  t  ; 


( A I  1  —  10’) 


\ 


We  have 


[L]  lyl  =  lb] 


(All-11) 


Equation  (AII-11) ,  which  can  be  written  explicitly  as  follows 


112yl  +  122y2 


1  ,y  +1  0y„  +. . .+1  y  =  b 
nl  1  n2  2  an  n  n 


(AII-I2) 


can  be  solved  by  the  trivial  forward  substitution  method,  which  merely 
solves  the  set  of  equations  (All-12)  sequentially  by  direct  substitution 
of  the  previously  calculated  y.’s. 

Next  we  write  Lquation  (All-10)  explicitly  as 


U  X  +  u.  ,x  , 


•+U,  X  =  y. 
Ln  n  1 


.+u,  >:  =  v , 

2n  *  i  1 


it  is  now  obvious  that.  the  unknown  x.*s  can  be  so  Lvov!  sequent  i  a  i  i  \  bv 

L 

a  bacK-forward  substitution  method  starting  Irom  the  last  equation  inco 
/  T  in  Via  been  solved  in  1. qua  Lion  (A  11-12). 
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